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Abstract 

We derive the effective coarse-grained field equation for the scalar condensate 
of the linear sigma model in a simple and straightforward manner using linear 
response theory. The dissipative coefhcient is calculated at tree level on the basis 
of the physical processes of sigma-meson decay and of thermal sigma-mesons and 
pions knocking sigma-mesons out of the condensate. The field equation is solved 
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for hot matter undergoing either one or three dimensional expansion and cooling 
in the aftermath of a high energy nuclear collision. The results show that the time 
constant for returning the scalar condensate to thermal equilibrium is of order 2 
fm/c. 

1 Introduction 

Scalar condensates often appear in quantum field theories when a symmetry is sponta- 
neously broken. Prominent examples include the Higgs condensate and the chiral con- 
densate. The equilibrium behavior of these condensates as a function of temperature and 
density has been extensively studied in the context of cosmology and heavy ion colli- 
sions. However, much remains to be learned about how these condensates really evolve 
in out-of-equilibrium systems. 

In this paper, we study the dynamical evolution of the chiral condensate in the 0{N) 
linear sigma model and apply it to the expanding matter in a high energy heavy ion 
collision. We recall that the sigma field a represents the quark condensate qq in the 
sense that they both have the same transformation properties for = 4, corresponding 
to two flavors of massless quarks. At high temperatures, quarks and gluons exist in a 
deconfined, chirally symmetric phase. At some critical temperature of order 150 MeV a 
transition to a hadronic phase occurs. In this confined and symmetry-broken phase the 
quark condensate is nonzero. In a collision between large nuclei the beam energy must 
be high enough so that the matter reaches at least approximate thermal equilibrium at 
a temperature greater than the critical temperature. We derive and solve an effective 
coarse-grained equation of motion for the chiral condensate, or mean sigma field, starting 
at the critical temperature and continuing down to the temperature where the system 
loses its ability to maintain local equilibrium and freezes out. 

The equation of motion approach we espouse here was used by Linde in his pioneer- 
ing work on phase transitions in relativistic quantum field theory although he did 
not use this technique to analyze dissipative processes or fields out of equilibrium. Our 
derivation of a coarse-grained field equation is based on linear response theory [Q, where 
the connection between finite temperature averages and the evolution of an observable in 
real physical time is clear and straightforward. There is less clarity, in our opinion, with 
the Influence Functional ^ and closely related Closed Time Path 0] methods as used 
by many others in this context ||^, |^, 0, ^ ^ . Full nonlinear response theory is certainly 
equivalent to both of those methods, as indeed it must be since they all describe the same 
physics. However, in most cases practical calculations can only be performed when the 
deviation from thermal equilibrium is small. This is true of the Influence Functional and 
Closed Time Path methods as well as response theory. Still we are able to make several 
improvements to the most recent treatment of the linear sigma model. First, the com- 
plications associated with doubling of the field variables do not arise, nor those arising 
from paths in the complex time plane. Second, we make a direct connection with physical 
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response functions. It is clear that these response functions should be evaluated in the 
fully interacting equilibrium ensemble. Here we shall estimate the relevant response func- 
tions by including the physical processes of sigma formation and decay via two pions, and 
the scattering of thermal pions and sigma mesons from sigma mesons in the condensate. 
To our knowledge, the latter processes have not been studied before in this or any related 
context. The amplitudes for these processes are computed in tree approximation. 

In section 2 we derive the fundamental equation. In section 3 we calculate decay and 
scattering contributions to the dissipative term in the coarse-grained field equation for 
the condensate. In section 4 we solve the equation of motion in the hadronic phase of a 
high energy heavy ion collision. Conclusions, and extensions of this paper, are discussed 
in section 5. 

2 The Fundamental Equation 

The Lagrangian of the linear 0{N) sigma model involves the sigma field and a vector of 
(A^ — 1) pion fields, tt: 

£ = i {d,af + i {d.Tvf - \X {a' + tt^ - f^f , (1) 

where A is a positive coupling constant and /tt is the pion decay constant. In the vacuum 

the symmetry is spontaneously broken: the sigma field acquires a vacuum expectation 
value of (0| a\0) = /t^, excitations of the sigma field have mass = 2A/^, and the pion is a 
Goldstone boson since we neglect explicit chiral symmetry breaking here. The symmetry is 
restored by a second order phase transition at the critical temperature = 12/^/ (A^+2). 

When one is interested in what happens at a fixed temperature at and below Tc the 
sigma field is usually expressed as 

a{x,t) =v + a'{yi,t) , (2) 

where the thermal average of the sigma field at temperature T is {(j)eq — v so that 
{<^')eq — 0. The equation of motion for the fluctuating component a' is 

a' - VV = A/^ (v + a')-X{v + af -X{v + a') tt" . (3) 

We now allow the sigma field to be slightly out of equilibrium, and write instead: 

a{yi,t)^v + as{x,t)+(jf{x,t). (4) 

Here (a) = v + ag, where v denotes the equilibrium value as before, but the deviation 
has been split into two pieces: a slow part ag, whose average is nonzero, and a fast part 
CT/, whose average is zero. (Primes have been dropped for clarity of presentation.) The 
notation (a) denotes averaging over the space-time volume chosen for coarse-graining and 
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thus the precise division between the fast and slow parts of the field depends on the choice 
that is made. For example, one may include in the slow part only those Fourier compo- 
nents with wavenumber below some cutoff value |]TU|, [T^], or one may average the field 
fluctuations over some coarse-graining time |jl2[. We shall not delve into any details here, 
but only remark that if the results depend strongly on the coarse-graining technique then 
the procedure is not very useful in the given context. In any case, the slow part represents 
occupation of the low momentum modes by a large number of particles and so may be 
treated as a slowly varying classical field. The fast part is a fully quantum field. The 
equilibrium ensemble average of an operator O is denoted by {0)eq and is characterized 
by the temperature T and the full Hamiltonian determined from the linear sigma model 
Lagrangian in the absence of as- When this equilibrium ensemble is perturbed by the 
presence of the field cXs the Hamiltonian is modified and the resulting (nonequilibrium) 
ensemble average is denoted by (O). 

It should be noted at this point that we have not allowed for a nonzero ensemble 
average of the pion field. Such a nonzero average is usually referred to as a disoriented 
chiral condensate, or DCC One could certainly allow for a DCC and carry through the 
following computations in a straightforward way. However in thermal equilibrium a DCC 
never develops (although it can arise from thermal fluctuations in a small system, but 
even then the probability is small [13||). This is in contrast to the scalar condensate, whose 
value is zero above the critical temperature and becomes nonzero below it. (A very small 
pion mass, too small even to affect the equation of state or correlation functions at the 
temperatures of interest, will still tilt the system towards a unique vacuum.) Thus we 
choose to focus on the behavior of the a field. 

Let us average the sigma field equation over time and length scales large compared 
to the scales characterizing the quantum fluctuations of the fields aj and tt but short 
compared to the scales typifying as- It is an assumption that such a separation exists, 
but in any given situation it can be verified or refuted a posteriori. Since (a/) = we 
obtain 



as - V' 



as 



\{v + as) f^-iv + asY 



3{aj) 



A(aJ) 



(5) 



The full ensemble averages are 



6{a') 



{a}),^ + 5{a}) , 



(6) 
(7) 
(8) 



where the deviations are caused by as and are generally proportional to as to some 
positive power. Equation (5) must be satisfied even when as vanishes. That determines 
the equilibrium value of v to be 



if T>Tr 



(9) 
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or 



J TT 



f/eq 



{7r^),,-{aj),Jv if r<T, 



(10) 



We are interested in the second of these because it represents the low temperature 
symmetry-broken phase. To first approximation in either (i) a perturbative expansion 
in A or (ii) an expansion in 1/A^ [Q, which are the usual approximations for the sigma 
model, the field fluctuations are 



f/eq 



d^p 1 
(2^£^ 



(11) 



where E^j = ^Jml + p^ and nsiE/T) = 1/ {exp{E/T) — 1) is the Bose-Einstein distribu- 
tion. Note that the sigma mass is still to be determined. The pion fluctuations are 



{^')eq 



{N-1 



'l2 



;i2) 



The latter follows from the fact that there are — 1 Goldstone bosons. Corrections to 
these formulae come from interactions not included in the effective mass. The sigma mass 
is very large at zero temperature and vanishes at Tc. The term {aj)eq/v is not zero on 
account of the cubic self-coupling of the a with coupling coefficient Xv. It does have a 
finite limit as f — *■ 0. However, once the approximations (11-12) have been made it is not 
legitimate to keep this term because it is one higher power in A and/or 1/A^, and keeping 
it would violate the 0{N) symmetry. Therefore the limits are: 



T 

T - 







/2 - (iV - 1)TV12 
f^-{N + 2)TV12 



(13) 



where the formula for was given earlier. It is, of course, consistent with the above 
expression. 

Using Eqs. (6-8) and (10) in (5), the coarse-grained field equation for a.^ can now be 
written as 



(To 



-A 



3val + + s{aj) + {v + a,)(35(aj) + ^(Tr^)) 



(14) 



Note that here the system is assumed to be at a fixed temperature so that v{T) is a 
constant. If the temperature is allowed to vary then time derivatives of v must be included; 
see the next section. Equation (14) shows that the sigma mass and the equilibrium value 
of the scalar condensate are related by = 2Af ^. They are temperature dependent and 
determined self-consistently from the formula 



2A 



d^p 1 
(2^^ 



(15) 
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The limits T — and T ^ Tc are readily obtained from those of v. An interpolating 
formula which connects these two limits is: 



V 



2A/2 /2 1 - [3/(iV + 2)] (TVT2) (1 - TVT2) 



(16) 



This is useful when studying solutions of the equation of motion numerically. 

The deviations in the fluctuations due to the presence of give rise to a renormaliza- 
tion of the parameters in the equation of motion of ag but, more importantly, they lead to 
dissipation. Energy can be transferred between the field cr^ and the fields cx/ and tt. These 

[TH 



as 



deviations in the fluctuations may be computed using linear response theory 
long as as is small. Here small means in comparison to either v (most relevant at low tem- 
perature) or to \J {(yf) (most relevant at high temperature). Technically, linear response 
theory is an expansion in powers of the Hamiltonian coupling the out-of-equilibrium field 
(Ts to the other modes of the system, and guarantees that symmetries of the theory are 
respected. This piece of the Hamiltonian consists of positive powers of as and so becomes 
smaller and smaller with decreasing departures from equilibrium. The coupling between 
the slow modes and the fast modes are determined straightforwardly from the potential 
to be 



Hs 



X (vas 



1^2 
2^s 



TT 



(17) 



and 



X 



al + + 3t;2 - fl + |(o^s + 2f )(T/ + a 
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In order to apply linear response analysis we need some initial conditions. In a nuclear 
collision, or in the early universe for that matter, it is assumed that the system reaches a 
state of thermal equilibrium at some negative time, that the system expands and cools, 
and at time t = the system is at the critical temperature. This implies the initial 
condition crs(x, t = 0) = 0. (In a finite volume one may wish to consider an ensemble 
of initial values chosen from a canonical distribution [|13].) From the standard theory of 
linear response one immediately deduces that 



eq 



6{a]{x)) = iX dt' J d'x'[{3v^-f^)as{x') + 3va^^{x') + a^^{x')){[af{x'),a]{x) 
+ i3X f dt' j dV (vas{x') + \al{x')) i\^a){x'),a']{x) 





+ zA^ dt' j d^x'as{x') 



aj{x'),a]{x) 



eq 



and 



6{n\x)) =iX I dt' j dV (vas{x') + \al{x')) {[tz\x'),'k^{x)\) 



'eq ■ 



(19) 



(20) 



The response functions are just the commutators of powers of the field operators evaluated 
at two different space-time points in the unperturbed (equilibrium) ensemble. We should 
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emphasize that the unperturbed ensemble does include all interactions among the fast 
modes and makes no approximation regarding the strength of these interactions. Insertion 
of these deviations into Eq. (14) represents the fundamental equation of this paper. For 
very small departures from equilibrium it is sufficient to keep only those terms which are 
linear in ag- The high temperature, symmetric phase is obtained by setting v = 0. 

The time-delayed response of the fast modes to the slow one evident in eqs. (19-20) 
has two effects: The sigma mass and self-interactions are modified, and dissipation occurs. 
These effects may be seen by expanding the slow field (Ts(x',t') in a Taylor series about 
the point (x, t) in Eqs. (p!9H20|) (although such an expansion is not required) . Terms with 
no derivative or an even number of derivatives either renormalize existing terms in the 
equation of motion or add new nondissipative ones, such as agV^Cg and V^V^cr^. Terms 
with an odd number of derivatives are explicitly dissipative. Examples are ag, (Ts&s, and 
V^a, PI. 

Perhaps the closest analysis to ours is due to Rischke [|^. The differences may be 
summarized thusly: First, Rischke used the infiuence functional method which is 
closely related to the closed-time-path method [Q , for deriving the equations of motion of 
the classical field. We use linear response theory. Since the former techniques ultimately 
rely on a perturbative expansion in terms of as anyway, one might as well employ linear 
response theory to begin with. Linear response theory is quicker, easier to use, and more 
intuitive. Second, we write down the field equation for as, which is the deviation of 
the scalar condensate from its equilibrium value v. Rischke writes down the equation 
of motion for cr = w -|- as- Of course these approaches are equivalent. Third, Rischke's 
response functions are valid for free fields only. For example, a response function from 
Eqs. (p!9|-pOD has the form ([0^(x', t'), 0^(x, i(;)])eg. For free fields this may be written 
2[D'^{x' — x) — D'^{x' — x)] where the D are Wightmann functions/propagators and the 
subscript indicates whether t' is greater than or less than t. Using our approach it is clear 
that these response functions should be evaluated in the fully interacting ensemble (but 
unperturbed by as). Taken together with the second difference this constitutes a major 
improvement over Rischke's analysis. This will be discussed in more detail in the next 
section. Fourth, Rischke employed a particular coarse graining technique by separating 
soft and hard modes according to the magnitude of the momentum. We have left the 
coarse graining technique open. Finally, with a view towards the formation of disoriented 
chiral condensates, or DCC, Rischke allowed for slow classical components of the pion 
field. We have not done so here, mainly to keep the analysis as simple and direct as 
possible, although it could easily worked out in the same way as the sigma field. 

3 Estimating the Response Functions 

One's first inclination might be to evaluate the response functions using free fields. Sup- 
pose, for example, that as is so slowly varying in space and time that to first approximation 
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it can be taken outside the integration over y and t' . Then it is a simple exercise to show 
that for a free field with mass m and energy E one gets 



f^dt' j dV([02(x',O,0'(x,t)])^^ = z j^*rft'/ ^-^(2nB(i?/T) + l)sin[2E(t-t')] • 

(21) 

The temperature-independent piece is a vacuum contribution and may be dropped for the 
present discussion. In the case that m = the momentum integral is done trivially with 
the following result: 

27rT 



47r2 



t 

ds 







tanh(27rTs) 

In the limit that t becomes large compared to 1/2ttT this approaches the asymptotic 
value iTt/2'ir. This corresponds to the first term of a Taylor expansion of as and further 
terms in the series bring in powers of {t — t') yielding dissipative coefficients that grow as 
t^, and so on. This is clearly unacceptable. The origin of this problem can be traced 
to the inadequacy of evaluating the response functions in the free field limit. Indeed, 
these response functions are closely related to the shear and bulk viscosities via the Kubo 
formulas which express those quantities in terms of ensemble averages of commutators of 
the energy-momentum tensor density operator at two different space-time points. It is 
known that determination of the viscosities requires summation of all ladder diagrams at 



a minimum [p!q] . 



An alternative is to make a harmonic approximation [0 instead of a Taylor series 
expansion which, in this situation, can be formulated as: 

Csit — s) ~ o"s(t) cos{m„s) — d-s(t) sm{m„s)/ma- . 

For the applications we have in mind in this paper the deviation of the scalar condensate 
from its equilibrium value is not expected to oscillate significantly; rather, it is expected 
to decay exponentially to zero. Therefore the harmonic approximation does not solve the 
problem. 

Rather than attempt a sophisticated evaluation of the response functions in this paper 
we shall proceed to estimate them based on direct physical reasoning. There are two 
obvious mechanisms for adding or removing sigma mesons from the condensate: (1) decay 
into two pions or the reverse process of formation via two pion annihilation, and (2) a 
meson from the thermal bath elastically scattering off a sigma meson and knocking it out 
of the condensate. The first of these processes is included in all analyses of the linear 
sigma model; the second of these processes has not been studied in the literature to our 
knowledge. We shall evaluate them at the tree level. This means that the parameters in 
the Lagrangian are to be fitted to experimental data at the tree level for consistency. Let 
us examine each of these in turn. 
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The decay rate of a sigma meson into two pions in the sigma's rest frame is |]T7 

(Ar-l)AV 



Svr TJlrr IQlT 



-Am, , (22) 



when account is taken of the relation = 2Af^. The decay rate for a sigma meson at 
rest in the finite temperature system is Bose-enhanced by a factor of [1 + nB{mcr/2T)]'^, 
where is the Bose-Einstein occupation number with the indicated argument of its 
exponential. The rate for two pions to form a sigma meson at rest is obtained from 
detailed balance by multiplying Fo-^TTTr by n^(mo-/2T). Therefore the net rate is 

r^TTTT = —rr: Am<^ [1 + ub) - = — — Xm^[l + 2nB\ 

iovr L J iOTT 

= ^^~^^ Am,coth(mj4T). (23) 

(The argument of hb is the same as above.) This leads to the dissipative term T^t^t^&s in 
the equation of motion . It agrees with eq. (80) of Rischke |g] . 

Scattering of a thermal boson 6 off a sigma meson with negligibly small momentum 
and which is considered to be a component of the background field as can be studied by 



evaluating the self-energy |19 



r dEpnB{E/T)^sfr\s). (24) 

J rrih 



7T J rrif, 

Here E and p are the energy and momentum of the boson b, s = + ml + 2m„E, and 
fcrij is the forward scattering amplitude. The normalization of the amplitude corresponds 
to the standard form of the optical theorem 

All 

_lni/(=-)(s), (25) 

'?cm 

where gem is the momentum in the cm frame. Because the cross section is invariant under 
longitudinal boosts the scattering amplitude transforms as follows: 

The imaginary part of the self-energy 

Imn,fe = --^/ dEp''nBiE/T)a^,is) (27) 

determines the rate at which the field decays |[18|| : 

r^fe = -ImU^b/m^ . (28) 
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The applicability of this expression is limited to those cases where interference between 
sequential scatterings is negligible. 

First consider pion scattering. We will calculate to tree level only and suppose that 
the parameters of the theory are adjusted to reproduce low energy experimental data at 
this same level. The invariant amplitude is 



M 



-2A 



1 



1 

u 



(29) 



The s, t and u are standard Mandelstam variables satisfying s + t + u = 2m^. Note that 
the forward scattering amplitude evaluated at threshold, A1(s = m^,t = 0), vanishes 



in accordance with Adler's Theorem [^. The differential cross section in the center-of- 
momentum frame is 

l-^l' (30) 



The total cross section is given by 



1? 



a 



+ 



is' + If 
2(s'2 - 3s' - 1 



+ 



6471^3 ■ 



9s' 



{S' - 1)3 



S 

In [s'{2 



/2 



6(s 



12 



is' - ly 



In 



J2 



s' + l 



s' 



(31) 



where s has been scaled by the sigma mass: s' = s/m^. The total cross section has a 
branch point and a pole at s = 2m^ due to the u channel exchange of a sigma meson. 
Going beyond the tree level is necessary to incorporate the finite width of the sigma meson 
and regulate the singularity. Indeed, the tree approximation is not reliable at high energy. 
For example, it is well known that unitarity is violated in tttt scattering at energies of order 
of one to two times m„ [^. Therefore we are only allowed to use the low energy limit, 
which is quite acceptable when T <^ rria- In this limit 



cr 



(32) 



where p is the pion momentum in the sigma rest frame. Notice the workings of Adler's 
Theorem here: According to that theorem the forward scattering amplitude / evaluated 
in the rest frame of any target particle must vanish like as p 0. Notice also that 
the total cross section cannot be obtained from the imaginary part of / because we are 
essentially using a Born approximation. In the limit that T Tc we have the opposite 
situation where rria -C T. Then the three-point vertices do not contribute {Ai = — 2A in 
place of (29)) and so use of 

(T = — (33) 



Atts 



may be considered more appropriate. 
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The contribution to the imaginary part of the self-energy is readily determined in 
these two Umits. At low energy/temperature 

imn.. = -:^:^C(7)A2^, (34) 

and at high energy /temperature 

ton„ = -^. (35) 

These are the contributions from a single pion and must be multiplied by iV — 1 to factor 
in all pions. For the accuracy required in this paper we can construct a Fade approximant 
to represent these results. 

Notice that this contribution to F diverges as T Tc on account of division by m„. 
This ought to come as no surprise since one is approaching a critical point where certain 
fluctuations diverge. 

The above analysis may be repeated for a sigma meson knocking another out of the 
condensate. The invariant amplitude for aa — > aa is 



M = -6A 



1 + 2,ml \ — ^-^ + '-^ + 



s — u — rri^ t — ^ 



(37) 



reflecting the symmetry in the s, t, u channels. The total cross section is 

Stts + 18 I2{s'^ -?,s'-l) , , , 

= 3 7 + 3 n - 73 7773 ^773 77 " 3) . (38) 



9A2 \s'-l) s'-2> (s'-4)(s'-2)(s'-l) 
This expression has no singularities because s' > 4. The low energy limit is 

9 

'^^f^, (39) 
2n mi 



which gives rise to the low temperature limit 

lmn^<, = --^A'T'e-"''/^ . (40) 

In the high energy limit only the four-point vertex contributes, giving 

9 A^ 
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which gives rise to the high temperature hmit 

3 



Imll^, = ---X'T\ (42) 
o47r 



This expression is identical to that calculated by Jeon and Rischke in a symmetric 
0^ model to which it can be compared. The two limits can be combined in a Pade 
approximant as 

9 A^T^ 

Imn.. ~ [96 + (e-^/^ - 1)] ' ^^^^ 

which is useful for numerical computations. 

The total rate is obtained by addition of all components, namely: 

= -ImU^^/ma , (45) 



where 
and 



= -ImU^^/m^ . (46) 



Other works, such as [^, do not include the latter two scattering contributions arguing 
that they are of order A^, while Fo-tttt is of order A. While this is correct at low T when 
rricr is large, it becomes questionable near to the critical temperature where mo- is small 
(see the discussion of Fig. 1 below). 

It should be noted that the scattering contributions diverge at the critical temperature 
on account of division by the vanishing sigma mass. This may be a signal of the breakdown 
of the use of tree-level scattering amplitudes and requires further investigation. 



4 Solution in an Expanding Fireball 

In this section we shall study solutions to the coarse-grained field equation in several 
limits. First, suppose that the volume and temperature are fixed in time but that the 
system is slightly out of equilibrium in the sense that ag 7^ 0. The field equation is then 

as + T{T)crs + ml{T)as = -A {3v{T)al + a^) . (47) 

where F is the sum of decay and scattering terms as given in the previous section. When 
|(Js| <^ V this equation can be linearized. It is equivalent to a simple damped harmonic 
oscillator. The system is underdamped if > F/2 and overdamped if < F/2. 
We choose A = 18 so that the sigma mass in vacuum is 6/,^, corresponding to the s- 
wave resonance observed in the vrvr channel. In figure 1 we plot rricr and the individual 
contributions to F as functions of temperature T. The field is overdamped when T > O.STc. 



12 



Next consider the expansion of the system created in a high energy nuclear colhsion. 
If the beam energy is high enough it will form a quark-gluon plasma with temperature 
greater than Tc. This "fireball" will expand and cool, eventually reaching T^. At this 
moment, say at time t = t^, the initial conditions for the coarse-grained field must be 
specified. We will assume, for sake of illustration, that the system is locally uniform so 
that spatial gradients are unimportant. Depending on whether the system is expanding 
spherically or only longitudinally along the beam axis (Z) = 3 or 1, respectively) we obtain 
the modified equation of motion 

as+v+ ji&s + v) + T{T)a, + ml{T)a, = -A {3v{T)a^^ + a^) . (48) 

In this situation t is really the local, or proper, time and the term proportional to D/t 
may be thought to arise from either the d'Alembertian or from a volume dilution term 
\V{t) /V{t)\as [^, ^ analogous to the Hubble expansion ||2^. Perhaps the easiest way to 



s 



obtain this equation is to derive the equation of motion for the total condensate a = v + a. 
and then make the substitution. Most authors actually solve the equation of motion for 
CT, but this is a matter of taste. Note that v{T{t)) is the instantaneous value of the 
equilibrium condensate and that is why no potential for it appears above. 

In addition to the coarse-grained field equation one needs to know how the local 
temperature evolves with time. It must be assumed that it changes slowly enough so that 
local equilibrium can be maintained by the rapidly fluctuating pion and sigma fields. In 
principle one ought to solve the equation 

de/dt = -Dw/t (49) 

where e is the energy density, w = e+P is the enthalpy and P is the pressure. Rather than 
working out a detailed description of the equation of state, which is really dominated by 
the rapidly fluctuating thermal fields, we simply assume a free massless gas of pions where 
the pressure is proportional to T^. The sigma meson has a mass small compared to T 
only very near the critical temperature, and its effect is therefore generally unimportant. 
As a consequence, the temperature falls with time according to the law 

m = j . (50) 

A reasonable numerical value for tc is 3 fm/c |^6|, |2^, The back reaction of as 



on the time evolution of the temperature is thereby neglected. This is a quite reasonable 
approximation because very little of the total energy resides in the field as- 

The equation of motion may be solved by numerically evolving an analytic solution 
in the neighborhood of tc- As t ^ tc the equilibrium condensate f — 0, however Eqs. 
(16) and (50) indicate that the v and i) are singular, as are the scattering contributions 
to the width F. The analytic behavior of cTg as t tc is uniquely determined by requiring 
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that the derivatives of ag exactly cancel these singularities, while ag 0. The result is 
displayed in figure 2 for D = 1 and in figure 3 for D = 3. The deviation from equilibrium 
as is maximal less than 1/2 fm/c after the critical temperature is passed, and it dies away 
with a time scale of about 2 fm/c. There is a hint of oscillatory motion in the solutions, 
but basically they are overdamped. 

5 Conclusion 

In this paper we have studied the dynamical evolution of the scalar condensate in the 0{N) 
linear sigma model in out-of-equilibrium situations. Our method is based on the equation 
of motion. Dissipation arises because of the response of the correlation functions of the 
fast modes to the slow modes of the fields. This is treated with standard linear response 
theory. These response functions should be computed exactly and used in the resulting 
dissipative, coarse-grained equation of motion. However, such explicit computations are 
generally not possible to do. Therefore, we identified the physical mechanisms responsible 
for the dissipation and estimated the corresponding response functions based on them. 
These mechanisms include the decay of sigma mesons in the condensate, and the knockout 
of sigma mesons in the condensate due to scattering with thermal sigma mesons and pions. 
To our knowledge, the latter physical mechanisms have not been studied before. 

We then studied the dynamical evolution of the condensate in heavy ion collisions, 
after the phase transition from quark-gluon plasma to hadrons, and allowing for either one 
or three dimensional expansion of the hot matter. These showed that thermal equilibrium 
was reestablished rather rapidly, with a time constant of order 2 fm/c. Clearly, much more 
could be studied along these same lines, including the formation and fate of disoriented 
chiral condensates (DCC). 

The method we used in this paper is very general, and may be applied to other 
theories, including nuclear matter, QCD, and electroweak theory. Such work is underway. 
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Figure 1: The sigma mass and the decay and scattering contributions to the width of the 
sigma meson as functions of temperature from to Tp. 
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Figure 2: The temporal evolution of the deviation of the scalar condensate from its equi- 
librium value V in units of MeV for a one dimensional expansion of hot matter produced 
in a high energy nuclear collision. 
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Figure 3: The temporal evolution of the deviation of the scalar condensate from its 
equilibrium value v in units of MeV for a three dimensional expansion of hot matter 
produced in a high energy nuclear collision. 
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